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Abstract 

We consider different aspects of the problem of cosmological singular- 
ity such as the BKL oscillatory approach to singularity, the new features 
of the cosmological dynamics in the neighbourhood of the singularity in 
multidimensional and superstring cosmological models and their connec- 
tions with such a modern branch of mathematics as infinite-dimensional 
Lie algebras. Besides, we consider some new types of cosmological singu- 
larities which were widely discussed during last decade after the discovery 
of the phenomenon of cosmic acceleration. 

1 Introduction 

Many years ago, in conversations with his students, Lev Davidovich Landau used 
to say that there were three problems most important for the theoretical physics: 
the problem of cosmological singularity, the problem of phase transitions and the 
problem of superconductivity [lj. Now we know that the great breakthrough was 
achieved in the explanation of phenomena of superconductivity 2 and phase 
transitions [3 . The problem of cosmological singularity was widely studied 
during the last 50 years and many important results were obtained, but it still 
conserves some intriguing aspects. Moreover some quite unexpecting facets of 
the problem of cosmological singularity were discovered. 

In our review, published 10 years ago [3] in the issue of this journal, dedicated 
to the 90th anniversary of L D Landau birth, we had discussed some questions 
connected with the problem of singularity in cosmology. In the present paper 
we shall dwell on relations between well-known old results of these studies and 
the new developments of this area. 

To begin with, let us remember that R Penrose and S Hawking [5] proved 
the impossibility of indefinite continuation of geodesies under certain conditions. 
This was interpreted as pointing to the existence of a sigularity in the general 
solution of the Einstein equations. These theorems, however, did not allow for 
finding the particular analytical structure of the singularity. The analytical be- 
haviour of the general solutions of the Einstein equations in the neighborhood 
of singularity was investigated in the papers by E M Lifshitz and I M Khalat- 
nikov 00 |S] and V A Belinsky, E M Lifshitz and I M Khalatnikov [9"l ITUI fPT] . 
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These papers revealed the enigmatic phenomenon of an oscillatory approach to 
the singularity which has become known as the Mixmaster Universe |12j . The 
model of the closed homegeneous but anysotropic universe with three degrees 
of freedom (Bianchi IX cosmological model) was used to demonstrate that the 
universe approaches the singularity in such a way that its contraction along two 
axes is accompanied by an expansion with respect to the third axis, and the axes 
change their roles according to a rather complicated law which reveals chaotical 
behavior QUI ED EH [H] . 

The study of the dynamics of the universe in the vicinity of the cosmological 
singularity has become a explodingly developing field of modern theoretical and 
mathematical physics. First of all we would like to mention the generalization 
of the study of the oscillatory approach to the cosmological singularity in mul- 
tidimensional cosmological models. It was noticed [T5] that the approach to the 
cosmological singularity in the multidimensional (Kaluza-Klein type) cosmologi- 
cal models has a chaotic character in the spacetimes whose dimensionality is not 
higher then ten, while in the spacetimes of higher dimensionalities a universe af- 
ter undergoing a finite number of oscillations enters into monotonic Kasner-type 
contracting regime. 

The development of cosmological studies based on superstring models has 
reavealed some new aspects of the dynamics in the vicinity of the singularity |17j . 
First, it was shown that in these models exist mechanisms of changing of Kasner 
epochs provoked not by the gravitational interactions but by the influence of 
other fields present in these theories. Second, it was proved that the cosmological 
models based on six main superstring models plus D = 11 supergravity model 
exhibit the chaotical oscillatory approach towards the singularity. Third, the 
connection between cosmological models, manifesting the oscillatory approach 
towards singularity and a special subclass of infinite-dimensional Lie algebras 
[IB] - the so called hyperbolic Kac-Moody algebras was discovered. 

Another confirmation of the importance of the problem of singularity in gen- 
eral relativity has come from the observational cosmology. At the end of 90th 
the study of the relation between the luminosity and redshift of supernovae of 
the type la has reavealed that the modern Universe is expanding with an accel- 
eration [19] . To provide such an acceleration it is necessary to have a particular 
substance which was dubbed "dark energy" [2D] . The main feature of this kind 
of matter is that it should have a negative pressure p such that p + 3p < 0, where 
p is the energy density. The simplest kind of this matter is the cosmological con- 
stant for which p = —p. The so called standard or ACDM cosmological model 
is based on the cosmological constant, whose energy density is responsablc for 
roughly 70 percents of the general energy density of the Universe, while the rest 
is occupied by a dust-like matter both baryonic (approximately 4 percents) and 
dark one. This model is in a good agreement with the observations, however 
other candidates for the role of dark energy are intensively studied and new ob- 
servations can give some surprises alredy in the nearest years. First of all let us 
notice that some observations [H] point out on the possible existence of the so 
called superacceleration, which is connected with the presence of the phantom 
dark energy [22j . characterized by the inequality p < — e. The universe, filled 
with such kind of dark energy, at some conditions can encounter a very partic- 
ular cosmological singularity - the Big Rip [23j . When an expanding universe 
encounters this singularity it has an infinite cosmological radius and an infinite 
value of the Hubble variable. Earlier the possibility of this type of singularity 
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was discussed in paper [24] . 

Study of different possible candidates for the role of dark energy has stim- 
ulated the elaboration of general theory of possible cosmological singularities 
[23 [Ml W\ US \M ■ Tt is remarkable that while the "traditional" Big Bang or 
Big Crunch singularities are associated with vanishing size of a universe, i.e. 
with a universe squeezed to a point, these new singularities occur at finite or 
infinite value of the cosmological radius. The physical processes taking place in 
the vicinity of such singularities can have rather exotic features and their study 
is of great interest. Thus, we see that the development of both the theoreti- 
cal and the observational branches of cosmology has confirmed the importance 
of the problem of singularity in general relativity mentioned by Landau many 
years ago. 

The structure of the paper is the following: in Sec. 2 we briefly discuss the 
Landau theorem about the singularity which was not published in a separate 
paper and was reported in the book [30] and in the review [B]; in Sec. 3 we 
shall recall the main features of the oscillatory approach to the singularity in 
relativistic cosmology; Sec. 4 will be devoted to the modern development of the 
BKL ideas and methods, including the dynamics in the presence of a massless 
scalar field, multidimensional cosmology, superstring cosmology and the corre- 
spondence between chaotic cosmological dynamics and hyperbolic Kac-Moody 
algebras; in Sec. 5 we describe some new types of cosmological singularities and 
in Sec. 6 we shall present some concluding remarks. 



2 Landau theorem about the singularity 

Let us consider the synchronous reference frame with the metric 

ds 2 = dt 2 -j a pdx a dxP, (1) 

where "/ a f3 is the spatial metric. L D Landau pointed out that the determinant g 
of the metric tensor in a synchronous reference system must go to zero at some 
finite time provided some simple conditions on the equation of state are satisfied. 
To prove this statement it is convenient to write down the — component of 
the Ricci tensor as 

I dKg _ I 
2 dt A' 

where K a p is the extrinsic curvature tensor defined as 



Rl = ~-^- L MK$, (2) 



K a0 = — (3) 

and the spatial indices are raised and lowered by the spatial metric j a p- The 
Einstein equation for i?g is 

fl§=l£-~T, (4) 

where the energy-momentum tensor is 

Tf = (p + p)uiU 3 - 5{p, (5) 

where p,p and Ui are the energy density, the pressure and the four- velocity 
respectively. The quantity Tg — \T in the right-hand side of equation ([4|) is 

T°-^T=±(p+3p) + (p+p)u a u a , (6) 
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and is positive provided 

p + 3p > 0. (7) 

Thus, from Eq. Q follows that 

Because of the algebraic inequality 

KiK«> l -{K2) 2 (9) 



we have 

dK° . 1 



dt 6 

or 



(K) 2 < o (io) 



aii 

dtKg ~ 6 v y 

If if" > at some moment of time, then if t decreases, the quantity -J^ decreases 
to zero within a finite time. Hence goes to +oo and that means because of 
the identity 

that the determinant g goes to zero (no faster than t 6 according to the inequality 
(fTTj) ). If at the initial moment <, then the same result will be obtained for 
the increasing time. The similar result was obtained by Raychaudhuri [31j for 
the case of dust-like matter and by Komar [32] . 

This result does not prove that it is inevitable that there exists a true physical 
singularity, which is one belonging to the spacetime itself and is not connected 
with the character of the chosen reference system. However, it played an im- 
portant role stimulating the discussion about the existence and generality of 
the singularities in cosmology. Notice that the condition of the energodomi- 
nance used for the proof of the Landau theorem appears also in the proof 
of the Penrose and Hawking singularity theorem [5] . Moreover, the breakdown 
of this condition is necessary for an explanation of the phenomenon of cosmic 
acceleration. 

The Landau theorem is deeply connected with the appearance of caustics 
studied by E M Lifshitz I M Khalatnikov and V V Sudakov [33J and discussed 
between them and L D Landau in 1961. Trying to construct geometrically 
the synchronous reference frame, one start from the three-dimensional Cauchy 
surface and design the family of geodesies orthogonal to this this surface. The 
length along these geodesies serves as the time measure. It is known that these 
geodesies intersect on some two-dimensional caustic surface. This geometry 
constructed for the empty space is valid also in the presence of dust- like matter 
[p = 0). Such matter, moving along the geodesies concentrates on caustics, 
but the growth of density cannot be unbounded because the arising pressure 
destroys the caustics 0- This question was studied by L P Grishchuk [34]. Later 



In an empty space the caustic is a mathematical, but not a physical singularity. This 
follows simply from the fact that we can always shift its location by changing the initial 
Cauchy surface. 
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V I Arnold, S F Shandarin and Ya B Zeldovich 35! have used the caustics for the 
explanation of the initial clustering of the dust, which though does not create 
physical singularities, is nevertheless responsible for the creation of the so called 
pancakes. These pancakes represent the initial stage of the development of the 
large scale structure of the universe. 



3 Oscillatory approach to the singularity in rel- 
ativist ic cosmology 

One of the first exact solutions found in the framework of general relativity 
was the Kasner solution [TB] for the Bianchi-I cosmological model representing 
gravitational field in an empty space with Euclidean metric depending on time 
according to the formula 

ds 2 = dt 2 _ t 2 P l dx 2 _ t 2 P2dy 2 _ t a Psd? 2 i ^3) 

where the exponents p\ , P2 and p^ satisfy the relations 

Pi+P2+P3=pI+pI+P 2 3 = 1- (14) 
Choosing the ordering of exponents as 

Pi <P2 < P3 (15) 

one can parameterize them as [1] 

-u 1 + u u(l + u) 

1 + u + u z 1 + u + u z 1 + u + u z 

As the parameter u varies in the range u > 1, p\,p2 and p% assume all their 
permissible values 

-g<J»i<0, 0<p2<|, |<ps<l. (17) 
The values u < 1 lead to the same range of values of ptiPiiPs since 

Pi (J^j =Pi(u), P2 Qj \ =P3(u), p 3 f^J =P2(u). (18) 

The parameter u introduced in early sixties has appeared very useful and its 
properties attract attention of researchers in different contests. For example, 
in the recent paper [35] a connection was established between the Lifshitz- 
Khalatnikov parameter u and the invariants, arising in the context of the Petrov's 
classification of the Einstein spaces [37] . 

In the case of Bianchi-VIII or Bianchi-IX cosmological models the Kasner 
regime (fT5 ]i . (TT¥j) ceased to be an exact solution of Einstein equations, however 
one can design the generalized Kasner solutions [7J-[Tl]. It is possible to con- 
struct some kind of perturbation theory where exact Kasner solution (Tl3]) . (fTi]) 
plays role of zero-order approximation while role of perturbations play those 
terms in Einstein equations which depend on spatial curvature tensors (appar- 
ently, such terms are absent in Bianchi-I cosmology). This theory of perturba- 
tions is effective in the vicinity of singularity or, in other terms, at t — > 0. The 
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remarkable feature of this perturbations consists in the fact that they imply the 
transition from the Kasner regime with one set of parameters to the Kasner 
regime with another one. 

The metric of the generalized Kasner solution in a synchronous reference 
system can be written in the form 

ds 2 = di 2 — (a 2 l a lp + b 2 m a mp + c 2 n a np)dx a dx^ , (19) 

where 

a = t Pl , b = t p ™, c = t p ™. (20) 

The three-dimensional vectors l,m,n define the directions along which the spa- 
tial distances vary with time according to the power laws (|20j) . Let pi — pi, p m — 
P2,Pn — P3 so that 

a~t p \ b~t p \ c-< P3 , (21) 

i.e. the Universe is contracting in directions given by vectors rh and h and is 
expanding along I. It was shown |10j that the perturbations caused by spatial 
curvature terms make the variables a, b and c to undergo transition to another 
Kasner regime characterized by the following formulae: 

a~ t p ', b~t p '\ c~ff3, (22) 

where 

, _ Iglj / _ 2 bll~P2 , _ P3-2|pi| 

Pl ~ l-2\ Pl \' Pm - l-2\ Pl \> Pn - l-2\ Pl \ ■ [Z6) 

Thus, the effect of the perturbation is to replace one "Kasner epoch" by 
another so that the negative power of t is transformed from the I to the rh 
direction. During the transition the function a(t) reaches a maximum and b(t) 
a minimum. Hence, the previously decreasing quantity b now increases, the 
quantity a decreases and c(i) remains a decreasing function. The previously 
increasing perturbation caused the transition from regime (I21|) to that (|22|) is 
damped and eventually vanishes. Then other perturbation begins grow which 
leads to a new replacement of one Kasner epoch by another, etc. 

We would like to emphasize that namely the fact that perturbation implies 
such a change of dynamics which extinguishes it gives us an opportunity to use 
perturbation theory so successfully. Let us note that the effect of changing of the 
Kasner regime exists already in the cosmological models more simple than those 
of Bianchi IX and Bianchi VIII. As a matter of fact in the Bianchi II universe 
there exists only one type of perturbations, connected with the spatial curvature 
and this perturbation makes one change of Kasner regime (one bounce). This 
fact was known to E M Lifshitz and I M Khalatnikov at the beginning of sixties 
and they have discussed this topic with L D Landau (just before the tragic 
accident) who has appreciated it highly. The results describing the dynamics 
of the Bianchi IX model were reported by I M Khalatnikov in his talk given 
in January 1968 in Henri Poincare Seminar, in Paris. John A Wheeler who 
was present there pointed out that the dynamics of the Bianchi IX universe 
represents a non-trivial example of the chaotic dynamical system. Later Kip 
Thorn has distributed the preprint with the text of this talk. 

Coming back to the rules governing the bouncing of the negative power of 
time from one direction to another one can show that they could be conveniently 
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expressed by means of the parameterization (116[) : 



Pi =Pi(u), Pm =P2(u), p n =P3(u) (24) 

and then 

Pi =pz(u - 1), p' m =pi(u- 1), p' n =ps(u- 1). (25) 

The greater of the two positive powers remains positive. 

The successive changes (l2"Sl . accompanied by a bouncing of the negative 
power between the directions / and m, continue as long as the integral part 
of u is not exhausted, i.e. until u becomes less that one. Then, according to 
Eq. (fT5|) the value u < 1 transforms into u > 1, at this moment either the 
exponent pi or p m is negative and p n becomes smaller one of the two positive 
numbers (p n = P2)- The next sequence of changes will bounce the negative 
power between the directions ft and I or ft and fh. Let us emphasize that 
the usefulness of the Landau-Khalatnikov parameter u is connected with the 
circumstance that it allows to encode rather complicated laws of transitions 
between different Kasner regimes (ffiJ)) in such simple rules as u — > u — 1 and 

Consequently, the evolution of our model towards a singular point consists of 
successive periods (called eras) in which distances along two axes oscillate and 
along the third axis decrease monotonically, the volume decreases according to 
a law which is near to ~ t. In the transition from one era to another, the 
axes along which the distances decrease monotonically are interchanged. The 
order in which the pairs of axes are interchanged and the order in which eras of 
different lengths follow each other acquire a stochastic character. 

To every (sth) era corresponds a decreasing sequence of values of the param- 
eter u. This sequence has the form Umax, Umax — 1, ■ ■ ■ , u m \ n , where u^ in < 1. 
Let us introduce the following notation: 

«£L = * W . = (26) 

i.e. k^ — [umL] (the square brackets denote the greatest integer < Umax)- 
The number defines the era length. For the next era we obtain 



a 



o+i) = _L 1.(8+1) 

max x (s) ' 



1 



As) 



(27) 



The ordering with respect to the length of k^ s ' of the successive eras (mea- 
sured by the number of Kasner epochs contained in them) acquires asymptoti- 
cally a stochastic character . The random nature of this process arises because 
of the rules (f2"6"l) - (f2"T)) which define the transitions from one era to another in 
the infinite sequence of values of u. If all this infinite sequence begins since 
some initial value Umax = k^ + , then the lengths of series fc^ ^ , fc^ 1 -* , . . . are 
numbers included into an expansion of a continuous fraction: 

+ = fc(°) + _ 1 , (28) 

We can describe statistically this sequence of eras if we consider instead of a 
given initial value Umax = k^ +x^ ^ a distribution of x ^ over the interval (0, 1) 
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governed by some probability law. Then we also obtain some distributions of 
the values of x^ s ' which terminate every sth series of numbers. It can be shown 
that with increasing s, these distributions tend to a stationary (independent 
of s) probability distribution w(x) in which the initial value x^ is completely 
"forgotten" : 

W ^ = (irW (29) 

It follows from Eq. (|29|) that the probability distribution of the lengths of series 
k is given by 

w < k > = ^WWy < 3 »> 

Moreover, one can calculate in an exact manner probability distributions for 
other parameters describing successive eras such as parameter S giving relation 
between the amplitudes of logarithms of functions a,b,c and logarithmic time 

Thus, we have seen from the results of statistical analysis of evolution in the 
neighbourhood of singularity [13j that the stochasticity and probability distri- 
butions of parameters arise already in classical general relativity. 

At the end of this section a historical remark is in order. The continuous 
fraction (T2g|) was shown in 1968 to I M Lifshitz (L D Landau has already left) 
and he immediately noticed that one can derive the formula for a stationary 
distribution of the value of x (f29|) . Later it becomes known that this formula 
was derived in nineteenth century by Karl F Gauss, who had not published it 
but had described it in a letter to one of colleagues. 



4 Oscillatory approach to the singularity: mod- 
ern development 

The oscillatory approach to the cosmological singularity described in the preced- 
ing section was developed for the case of an empty spacetime. It is not difficult 
to understand that if one considers the universe filled with a perfect fluid with 
the equation of state p = wp, where p is the pressure, p is the energy density 
and w < 1, then the presence of this matter cannot change the dynamics in the 
vicinity of the singularity. Indeed, using the energy conservation equation one 
can show that 

= Po = Po fall 
9 ~ (abc) w + 1 ~ fH-i ' { ' 

where po is a positive constant. Thus, the term representing the matter in 
the Einstein equations behaves like ~ \/t 1+w and at t — > is weaker than 
the terms of geometrical origin coming from the time derivatives of the metric, 
which behaves like 1/t 2 , let alone the perturbations due to the presence of 
spatial curvature, responsable for changes of a Kasner regime,which behave like 
l/i 2+4 l pi L However, the situation changes drastically, if the parameter w is 
equal to one, i.e. the pressure is equal to the energy density. Such kind of 
matter is called "stiff matter" and can be represented by a massless scalar field. 
In this case p ~ 1/t 2 and the contribution of matter is of the same order as 
main terms of geometrical origin. Hence, it is necessary to find a Kasner type 
solution, taking into account the presence of terms, connected with the presence 
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of the stiff matter (a massless scalar field) . Such study was carried out in paper 
[3"5] , It was shown that again the scale factors a, b and c can be represented as 
t 2pi , t 2p2 and t 2ps respectively, where the Kasner indices satisfy the relations: 



P1+P2+ P3 



1, pI+pI+pI 



9 

q , 



(32) 



where the number q 2 reflects the presence of the stiff matter and is bounded by 



q 2 < 



(33) 



One can see that if q 2 > 0, then exist combinations of the positive Kasner 
indices, satisfying the relations |3"2"|) . Moreover, if q 2 > | only sets of three 
positive Kasner indices can satisfy the relations ([3"2"]) . If a universe finds itself 
in a Kasner regime with three positive indices, the perturbative terms, existing 
due to the spatial curvatures are too week to change this Kasner regime, and 
thus, it becomes stable. That means, that in the presence of the stiff matter, 
the universe after a finite number of changes of Kasner regimes finds itself in 
a stable regime and oscillations stop. Thus, the massless scalar field plays 
"anti-chaotizing" role in the process of the cosmological evolution [38] . One can 
use the Lifshitz-Khalatnikov parameter also in this case. The Kasner indices 
satisfying the relations (|32p are conveniently represented as |38j 



Pi 

P2 
P3 



1 + u + u 2 ' 

1 + u 
1 + u + u 2 

1 + u 
1 + u + u 2 
2(1 +u + u 2 ) 2 

(U 2 - l) 2 



(1-/? 2 ) 1/2 ) 
1 + ^(1-(1-/3 2 ) 1/2 ) 



(34) 



The range of the parameter u now is — 1 < u < 1, while the admissible values 
of the parameter q at some given u are 



q 2 < 



(u 2 - l) 2 



2(1 + u + u 2 ) 



2\2 ' 



(35) 



One can easily show that after one bounce the value of the parameter q 2 changes 
according to the rule 



q' 2 =q 2 x 



1 



(l + 2pi)< 



(36) 



Thus, the value of the parameter q 2 grows and, hence, the probablity to find all 
the three Kasner indices to be positive increases. It confirms again the statement 
that after a finite number of bounces the universe in the presence of the massless 
scalar field finds itself in the Kasner regime with three positive indices and the 
oscillations stop. 

In the second half of eightees a series of papers was published [TS], where 
were studied the solutions of Einstein equations in the vicinity of singularity 
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for d + 1-dimcnsional spacetimes. The multidimensional analog of a Bianchi-I 
universe was considered, where the metric is a generalized Kasner metric: 



ds 2 =dt 2 -^t 2pi dx i2 , (37) 



i=l 

where the Kasner indices pi satisfy the conditions 

d d 



]>> = 5> 2 = 1. (38) 



i=l j=l 



In the presence of spatial curvature terms the transition from one Kasner epoch 
to another occurs and this transition is described by the following rule: the new 
Kasner exponents are equal to 



with 



Pi, ■■■,Pd = ordering of (qi,q 2 , ■ ■ .,%), (39) 



-Pi - P P2 



l + 2 Pl + P ' "* l + 2 Pl + P' 

Pd-2 _ 2pi +P +Pd-1 

l + 2 Pl +P' qd ^ ~ l + 2 Pl + P ' 
'i+P + Pd 
l + 2 Pl + P' 



qd-2 

2 Pl + P + Pd 

<?d = \ , „ —5- , (40) 



where 



d-2 
i=2 

However, such a transition from one Kasner epoch to another occurs if at least 
one of the numbers onjh is negative. These numbers are defined as 

a ijk = 2 Pi + Pu (i k, j^k). (42) 

For the spacetimes with d < 10 one of the factors a is always negative and, 
hence one change of Kasner regime is followed by another one, implying in 
such a way the oscillatory behaviour of the universe in the neighbourhood of 
the cosmological singularity. However, for the spacetimes with d > 10 exist 
such combinations of Kasner indices, satisfying Eq. (|38[) and for which all the 
numbers otijk are positive. If a universe enters into the Kasner regime with such 
indices, ( so called "Kasner stability region") its chaotical behaviour disappers 
and this Kasner regime conserves itself. Thus, the hypothesis was forward that 
in the spacetimes with d > 10 after a finite number of oscillations the universe 
under consideration finds itself in the Kasner stability region and the oscillating 
regime is replaced by the monotonic Kasner behaviour. 

The discovery of the fact that the chaotic character of the approach to the 
cosmological singularity disappears in the spacetimes with d > 10 was unex- 
pected and looked as an accidental result of a game between real numbers satis- 
fying the generalized Kasner relations (|40|) . Later it becomes clear that behind 
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this fact there is a deep mathematical structure, namely, the hyperbolic Kac- 
Moody algebras. Indeed, in the series of works by Damour, Henneaux, Nicolai 
and some other authors (see e.g. Refs. [T7]) on the cosmological dynamics in 
the models based on superstring theories, living in 10-dimcnsional spacetime 
and on the d + 1 = 11 supergravity model, it was shown that in the vicinity of 
the singularity these models reveal ocsillating bahaviour of the BKL type. The 
important new feature of the dynamics in these models is the role played by 
non-gravitational bosonic fields (p-forms) which are also responsable for transi- 
tions from one Kasner regime to another. For description of these transitions 
the Hamiltonian formalism p~2] becomes very convenient. In the framework of 
such formalism the configuration space of the Kasner parameters describing the 
dynamics of the universe could be treated as a billiard while the curvature terms 
in Einstein theory and also p-form's potentials in superstring theories play role 
of walls of these billiards. The transition from one Kasner epoch to another is 
the reflection from one of the walls. Thus, there is a correspondence between 
rather complicated dynamics of a universe in the vicinity of the cosmological 
singularity and the motion of an imaginary ball on the billiard table. 

However, there exists a more striking and unexpected correspondence be- 
tween the chaotical behaviour of the universe in the vicinity of the singularity 
and such an abstract mathematical object as the hyperbolic Kac-Moody alge- 
bras [17] . Let us explain briefly what does it mean. Every Lie algebra is defined 
by its generators hi, &i, fi, i = 1, . . . , r, where r is the rank of the Lie algebra, i.e. 
the maximal number of its generators hi which commutes each other (these gen- 
erators constitute the Cartan subalgebra). The commutation relations between 
generators are 

fj\ = &ijhi, 

[hi,ej\ = Aij = Aijej, 
[hi, fj] — —Aijfj, 

[h i ,h j ]=Q. (43) 

The coefficients Aij constitute the generalized Cartan r x r matrix such that 
An — 2, its off-diagonal elements are non-positive integers and Aij = for 
i ^ j implies Aji = 0. One can say that the are rising operators, similar to 
well-known operator L + = L x + iL y in the theory of angular momentum, while 
fi are lowering operators like L_ = L x — iL y . The generators hi of the Cartan 
subalgebra could be compared with the operator L z . The generators should 
also obey the Serre's relations 

(ad eiY'^ej = 0, 

(ad fif-^fj = 0, (44) 

where (a,dA)B = [A, B]. 

The Lie algebras Q (A) build on a symmetrizable Cartan matrix A have been 
classified according to properties of their eigenvalues: 
if A is positive definite, G(A) is a finite-dimensional Lie algebra; 
if A admits one null eigenvalue and the others are all strictly positive Q(A) is 
an Affine Kac-Moody algebra; 

if A admits one negative eigenvalue and all the others are strictly positive, G(A) 
is a Lorentz KM algebra. 



11 



There exists a correspondence between the structure of a Lie algebra and a 
certain system of vectors in the r-dimensional Euclidean space, which simplify 
essentially the task of classification of the Lie algebras. These vectors called 
roots represent the rising and lowering operators of the Lie algebra. The vectors 
corresponding to the generators and fa are called simple roots. The system 
of simple positive roots (i.e. the roots, corresponding to the rising generators 
ej) can be represented by nodes of their Dynkin diagrams, while the edges 
connecting (or non connecting) the nodes give an informartion about the angles 
between simple positive root vectors. 

An important subclass of Lorentz KM algebras can be defined as follows: A 
KM algebra such that the deletion of one node from its Dynkin diagram gives 
a sum of finite or affinc algebras is called an hyperbolic KM algebra. These 
algebras are all known. In particular, there exists no hyperbolic algebra with 
rank higher than 10. 

Let us recall some more definitions from the theory of Lie algebras. The 
reflections with respect to hyperplanes orthogonal to simple roots leave the 
systems of roots invariant. The corresponding finite-dimensional group is called 
Weyl group. Finally, the hyperplanes mentioned above divide the r-dimensional 
Euclidean space into regions called Weyl chambers. The Weyl group transform 
one Weyl chamber into another. 

Now, we can briefly formulate the results of approach [17 following the pa- 
per [3§] : the links between the billiards describing the evolution of the universe 
in the neigbourhood of singularity and its corresponding Kac-Moody algebra 
can be described as follows: 

the Kasner indices describing the "free" motion of the universe between the 
reflections from the wall correspond to the elements of the Cartan subalgebra 
of the KM algebra; 

the dominant walls, i.e. the terms in the equations of motion responsible for the 
transition from one Kasner epoch to another, correspond to the simple roots of 
the KM algebra; 

the group of refelections in the cosmological billiard is the Weyl group of the 
KM algebra; 

the billiard table can be identified with the Weyl chamber of the KM algebra. 

One can immagine two types of billiard tables: infinite such where the linear 
motion without collisions with walls is possible (non-chaotic regime) and those 
where reflections from walls are inevitable and the regime can be only chaotic. 
Remarkably, the Weyl chambers of the hyperbolic KM algebras are designed in 
such a way that infinite repeating collisons with walls occur. It was shown that 
all the theories with the oscillating approach to the singularity such as Einstein 
theory in dimensions d < 10 and superstring cosmological models correspond to 
hyperbolic KM algebras. 

The existence of links between the BKL approach to the singularities and the 
structure of some infinite-dimensional Lie algebras has inspired some authors to 
declare a new program of develpment of quantum gravity and cosmology [40] . 
They propose "to take seriously the idea that near the singularity (i.e. when 
the curvature gets larger than the Planck scale) the description of a spatial con- 
tinuum and space-time based (quantum) field theory breaks down, and should 
be replaced by a much more abstract Lie algebraic description" . 
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5 New types of cosmological singularities 



As was already mentioned in the Introduction, the development of the theoret- 
ical and observational cosmology and, in particular, the discovery of the cosmic 
acceleration has stimulated the elaboration of cosmological models where new 
types of singularities were described. In contrast to the "traditional" Big Bang 
and Big Crunch singularities, these singularities occur not at zero but at finite 
or even infinite values of the cosmological radius. The most famous of these 
singularities is, perhaps, the Big Rip singularity [231 HI] arising if the absolute 
value of the negative pressure p of the dark energy is higher than the energy 
density p. Indeed, let us consider a flat Friedmann universe with metric 

ds 2 - a 2 {t)dl 2 , (45) 
filled with the perfect fluid with the equation of state 

p = wp, w = const < — 1. (46) 
The dependence of the energy density p on the cosmological radius a is as usual 

C 



and the Friedmann equation in this case has the form 

a? C 

where C is a positive constant. Integrating Eq. (|48l) we obtain 



(47) 



(48) 




°(*) = K 2 + ^TTTT^ 1 ■ ( 49 ) 



It is easy to see that at the finite moment tu > to equal to 

3(1 + 10) 3(i+ m ) 
t R =t -^^-a * (50) 

the cosmological radius becomes infinite and the same occurs also with the 
Hubble variable - and, hence, with the scalar curvature. Thus, we encounter 
the new type of cosmological singularity, characterized by infinite values of the 
cosmological radius, its time derivative, the Hubble variable and the scalar cur- 
vature. Usually this singularity is called "Big Rip" singularity. Its properties 
have attracted an essential attention of researchers because of the fact that 
some observational data indicate that the actual value of the equation of state 
parameter w is indeed smaller than -1. 

There are also other types of cosmological singularities which can be encoun- 
tered at finite values of the cosmological radius (see e.g. [351 12S1 (2Z1 HE1 123 )• 
We shall consider here for illustration one type of singularities - the Big Brake 
singularity [25 . This singularity can be acchieved in a finite period of cosmic 
time and is characterized by a finite value of the cosmological radius, by a van- 
ishing first time derivative of the radius and by the second time derivative of 
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the cosmological radius tending to the minus infinity (an infinite deceleration) . 
Let us consider a perfect fluid with the equation of state 



P=~, (51) 
P 

where A is a positive constant. This fluid could be called "anti-Chaplygin" gas, 
because widely used the Chaplygin gas cosmological model [UJ is based on the 
equation of state p = —A/ p. The dependence of the energy density on the 
cosmological radius for the equation of state (j5Tj) is 



where B is a positive constant. When a is small p ~ 1/a 3 and behaves like dust. 
Then , when a — > a B 

a B = (~\ 6 (53) 

and the energy density tends to zero. The solution of the Friedmann equation 
in this limit gives 

o(t) = a B - C (t B - t)*, Co = 2- 7 / 3 3 5 / 3 (AB) 1 / 6 . (54) 

Now, one can easily check that when t — ► t B , 6 - * and a — > — oo. Thus, we 
indeed encounter the Big Brake cosmological singularity. 



6 Conclusions 

We have shown in this short review that the opinion expressed by L.D. Lan- 
dau many years ago concerning the importance of the problem of singularity 
in cosmology has proved to prophetic. The study of cosmological singularity 
has revealed the existence of the oscillatory behaviour of a universe when the 
curvature of spacetime grows, which in turn has a deep connection with quite 
new branches of modern mathematics. From the other hand the last successes 
of the observational cosmology has stimulated the development of various cos- 
mological models, which reveal new types of cosmological singularities, whose 
investigation from both physical and mathematical points of view can be very 
promising. 

This work was partially supported by the RFBR grant 05-02-17450 and by 
the grant LSS-1 157.2006.2. 
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